Thermodynamic properties of the d-density-wave order in cuprates 
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We solve a popular effective Hamiltonian of competing d-density wave and d-wave superconduc- 
tivity orders self-consistently at the mean-field level for a wide range of doping and temperatures. 
The theory predicts a temperature dependence of the d-density wave order parameter seemingly 
inconsistent with the neutron scattering and pSR experiments of the cuprates. We further calculate 
thermodynamic quantities, such as chemical potential, entropy and specific heat. Their distinct 
features can be used to test the existence of the d-density wave order in cuprates. 
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Unconventional charge and spin density wave orders 
were extensively investigated in correlated electron sys- 
tems 0. Recently, Chakravarty et al Q proposed that 
the pseudogap phenomena in the high T c superconduc- 
tors may originate from a hidden long-range order, d- 
density wave (DDW) (U ■ This state is also related to the 
staggered- flux state of Lee and Wen Gj , but the latter is 
dynamically fluctuating in their SU(2) gauge theory of 
the cuprates. According to Ref |^], the pseudogap is a 
consequence of the competition between two independent 
orders DDW and <i-wave superconductivity (DSC), which 
are transformable to each other in a 3 dimensional or- 
der parameter space and may coexist in the underdoped 
cuprates. A theory of DDW seems natural to account for 
a possible quantum critical point near the optimal doping 
level that marks the onset of the pseudogap, put forward 
by Tallon et al |5| by examining the data of photoemis- 
sion, thermodynamic and transport properties, etc. 

This DDW scenario has recently attracted much at- 
tention about its nature and experimental consequences 
H @i ||i H; @ ■ Some investigations || seem to indicate 
that various experiments in YBCO systems jll], [l2| |TJ] 
support this idea. The neutron scattering experiment 
jjr| shows that the elastic signal around the in-plane 
wave vector Q = (tt, it) in the underdoped YBCO ap- 
pears well above T c . The [j,SR experiment also con- 
firms that a small internal magnetic field appears above 
T c in the underdoped YBCO but below T c in the opti- 
mally doped samples. Since internal magnetic fields are 
very weak and spin fluctuations are too fast to couple 
with muon's spins, it is reasonable to attribute them to 
DDW. However, both experiments also show that the 
magnetic signal is enhanced when the temperature drops 
across T c . Such a behavior is not expected intuitively 
from the DDW picture, since the DSC and DDW orders 
compete each other. It thus becomes quite interesting 
to understand how this temperature-dependent puzzling 
behavior happens and, in particular, to see whether it 
can be understood in terms of the existing self-consistent 
DDW mean-field theory |, @, |, §, |l§. 

In this article, we shall study the temperature depen- 
dence of the DDW order at the mean-field level. We find 



that it cannot give the expected temperature dependence 
in the above experiments. Instead, our results show that 
the DDW order parameter is always suppressed when 
the temperature drops below T c , and that this is a ro- 
bust behavior, independent of the choice of parameters. 
Thermodynamic quantities (e.g. chemical potential, en- 
tropy, specific heat) are also investigated. Their distinct 
temperature dependences are discussed. 

Consider the following mean-field Hamiltonian |(| [?| ||, 

i0 

Hmf = {-Us - Vixtj)4 a c ja + h.c. 

- v 2 j2^(4A- c l c lO + h - c - 

ia (ij) (ij) 

where (ij) indicates summation over the nearest neigh- 
bors only. Aij and the imaginary part of Xij P ia Y the 
role of the DSC and DDW order parameters, respec- 
tively. They are related to the electron operators via 

Ay = (ct^ - ct ; c] T ) and Xij = (4* c j<r)- v i and 
V2 are positive in order to have a nonzero DDW and 
DSC order. We have used an effective hopping ampli- 
tude t c ff — tS with S the doping concentration and t the 
bare hopping amplitude to take account of the reduction 
of t near half-filling due to the strong Coulomb repulsion. 
Loosely speaking, the effective theory described above is 
equivalent to the fermion part of the slave-boson mean- 
field theory of the t- J model |L4|, |TJ] in which the J-term 
is decoupled into particle-hole and particle-particle chan- 
nels with different weights. 

The t- J model at half filling has a (hidden) local SU(2) 
symmetry |r|, which rotates (Imx, ReA, ImA) as a 3- 
vector. Thus the DDW(7r-flux) phase is degenerate with 
the DSC phase. Finite doping breaks this local SU(2) 
symmetry explicitly and favors DSC order 113 ], because 
the Fermi surface nesting is destroyed. Ubbens and Lee 
p5| showed that at finite temperatures the DDW (flux) 
state is stable only very close to half-filling. The bound- 
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FIG. 1: (A) Wddw vs temperature for various doping levels. 
From top to bottom, 5 = 0.05, 0.10, 0.14, 0.17, 0.18. (B) Adsc 
vs temperature for various doping levels. From bottom to top, 
8 = 0.10, 0.14, 0.18, 0.17. Curves of S = 0.18 are plotted with 
dashed lines for the aid of eyes. 



ary between the DDW(flux) and DSC phases is of the 
first order in nature. Thus there is no coexistence phase. 
In the model we are considering, V\ > V 2 is needed to 
have DSC and DDW coexist as pointed out in Ref. (7). 
Although the Heisenberg term equally favors DDW and 
DSC orders, the repulsion between the nearest sites fa- 
vors DDW over DSC. Hence the mean- field Hamiltonian 
above is reasonably postulated. 

Because the DDW order breaks the translational sym- 
metry, the Brillouin zone is reduced into one half and 
the operators (efef, Ck+Qf, c_ k ,, c_ k _Q^) are mixed to 
give two branches of Bogoliubov quasiparticle excita- 
tions: E(k)± = {(-/x ± W k f + (2V 2 Acb k ) 2 }^ 2 , where 

-{teg + ViRexhk and 
cos ky (eft is the tight 



W k = ^/el + (2V 1 lm X M 2 ,e k = 
4> k = cos k x — cos k y , 7ft = cos k x 
bond band energy). The corresponding self-consistent 
equations of Rex, Im% and A are: 



Rex 



1 

2N 



— tanh 



tanh 



pE k+ \ -fi + W k 



E k 



-fi -W k \ 

E k - y 



Imx 



2N^ 



tanh 



214 Imx 

wT 



tanh 



V 2 



-M + Wk 
E k + 



—fj, — Wk \ 

e, r 



A = 



S = 



k 

/E k _} x 2V 2 A4> : 



V 2 J E k+ V 2 



-I* - w k i 
E k - J 

where the summation is restricted within the reduced 
Brillouin zone and (3 = 1/T. Below we denote the energy 
gaps Adsc = 2F 2 A and Wddw = 2FiImx for the DSC 
and DDW orders, respectively. 

After solving the self-consistent equations at V\ = 0.38 
and V% = 0.25 with energy scale set as t = 1, we obtain 
the dependence of the Adsc and Wddw gap vs doping 8 
at zero temperature as shown in Fig. 1 of Ref. g]. Adsc 
begins to develop after S > 0.06 and reaches maximum 
at 5 w 0.18. Wddw also drops to zero around there. The 
phase diagram of temperature vs 6 is similar to Ref , 
and we shall not reproduce it here. 

There exists a coexisting region of both orders. How- 
ever, in this region the behavior Wddw vs temperature 
(T) is subtle (see Fig. [j]A). In comparison, Fig. |lB shows 
how Adsc varies with T. For very low doping ((5=0.05) 
where Adsc = 0, Wddw is monotonically enhanced 
when T is reduced. When the coexisting region is en- 
tered, Wddw becomes maximum around T c where Adsc 
starts developing. This feature is general for competing 
orders, which also happens in the competition of s and d- 
wave pairing orders |17|| . When either order develops, it 
generates a gap near the Fermi surface. Consequently, it 
becomes difficult for the other to form. When T slightly 
drops from T c , Adsc increases fast as (1 — T/Tc) 1 / 2 . 
Wddw loses more weight to DSC than it gains from low- 
ering temperature. When T drops well below T c , Adsc 
increases very slowly and Wddw changes little as well. 
In the underdoped region, Wddw still has a substantial 
residual value at T = K, which gets significantly re- 
duced and may even become zero near the optimal dop- 
ing. When 6 = 0.18 in Fig. [l]A, Wddw only exists in a 
small range around T c . In other words, something like 
the re-entry phenomenon occurs here with varying tem- 
peratures: Wddw vanishes at T = OK, begins to increase 
from a finite temperature to T c , and then decreases to 
zero again when T > T c . We never see Wddw develop- 
ing with decreasing T within the superconducting region. 
These are just opposite to what the neutron scattering 
and /iSR experiments indicated. 

We have to take it seriously why experimental signals 
are enhanced when T < T c . There seem to be only two 
possibilities. The first one is that these signals are really 
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FIG. 2: A. Chemical potential vs doping at zero tempera- 
ture. B. Chemical potential vs temperatures at characteristic 
doping levels. <5 = 0.05, 0.14, 0.20. 



related to the DDW order. Then a modified mean-field 
theory is needed for a mechanism that the two orders can 
somehow enhance each other around T c . Or these signals 
have other origins such as from spin, then they cannot 
be used as evidence for the existence of the DDW order. 

The behavior of chemical potential \i is also interest- 
ing. In Fig. |^A, we show the dependence of /i with S 
at zero temperature, fj, decreases with increasing S when 
5 < 0.06, increases slowly in the underdoped region and 
drops quickly in the overdoped region. This can be un- 
derstood as follows. The energy curve of quasiparticles is 
cone-shaped in the momentum space with a Fermi pocket 
near (n/2,n/2). After the onset of Adsc, Wddw drops 
faster with increasing 6, so the density of states (DOS) 
increases. This keeps \i roughly unchanged while increas- 
ing doping. In the underdoped region d/i/dn is small and 
becomes even negative, which means that the charge in- 
stability may develop here. Photoemission experiments 
p8[ show that /i is almost fixed at the undoped value 
upon increasing S in the underdoped region. Our result 
agrees with that qualitatively, but fj, is not fixed at the 
value of the undoped case. 

Fig. shows the temperature dependence of fi at 
three doping levels: low (nonsuperconducting) doping 
(5=0.05), underdoped(r5=0.14), and overdoped(r5=0.20). 
In the low doping region, fi increases with temperature 
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FIG. 3: Entropy per lattice site vs doping at fixed tempera- 
tures. From top to bottom, T=0.08, 0.04. 



increasing. Wddw is weakened by temperature while fi 
increases to fix the particle number. In the underdoped 
region, /i's behavior is subtle. It first drops when T < T c , 
then increases when T c < T < Tddw, and drops again 
when T > Tddw, where Tddw is the onset tempera- 
ture for DDW. This can be explained by the temperature 
dependence of Wddw (Fig. [y). Wddw is enhanced (or 
weakened) by increasing temperature within the range 
T < T c (or T c < T < Tddw)- Thus fi first drops and 
then increases. After T passes over Tddw, A* drops as 
an ordinary Fermi gas behaves. In the overdoped region, 
Wddw = 0. Thus \x is almost fixed when T < T c but 
drops when T > T c . 

The DDW order also has important effects on the en- 
tropy per site S vs S, as shown in Fig 0. The first curve 
is at the high temperature where Adsc = and only 
Wddw exists. In the underdoped region, S decreases 
when 5 is reduced, because Wddw reduces the low energy 
DOS. In the overdoped region, S drops when S increases, 
which is the standard Fermi liquid behavior. Thus S 
reaches maximum near the optimal doping. This agrees 
with experimental results of Loram et al . Very close 
to half filling, the hopping amplitude is reduced. As a 
result, the band width is reduced and DOS is enhanced. 
This effect tends to increase S. Simultaneously, the DDW 
order is enhanced by lowering doping, which has an ef- 
fect to decrease S. At high temperatures where Wddw is 
small, the first effect may overcome the second one and 
thus S increases when 5 decreases. This phenomenon 
is absent at the low temperature where Wddw is large, 
which is also shown in Ref |L9| . Let us increase S at the 
lower temperature as in the bottom curve. S increases 
at very low doping levels, since only Wddw exists. In 
the coexistence or underdoped region, S drops because 
Adsc develops. After passing the optimal point, Adsc 
decreases and thus makes S increase again. 

There is also specific heat anomaly at the onset of 
WddWi as shown in Fig 5. The jump of the specific 
heat coefficient y(T) = C(T)/T is large and can be com- 
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FIG. 4: Specific heat coefficient y(T) vs S in the underdoped 
and optimal region, from left to right S =0.08, 0.10, 0.17. 

pared with those at the onset of Adsc- The jump at the 
DDW transition is larger at the more underdoped side, 
while that at the DSC transition's behavior is just oppo- 
site. However, the former is not seen in experiments and 
is a difficulty for the DDW scenario, ft was argued that 
disorder removes the sharp transition and turns it into 
a smooth crossover in Ref. Ref. j^o) suggests that 
a negatively large fj, can weaken the jump by destroying 
the nesting of Fermi surface. However, from the Fig. 2B, 
\x increases rapidly when Tddw, and ImI is much 
smaller at Todw than at zero temperature, especially at 
the low doping region. It is still difficult to understand 
why [i can be negatively large. 

At last, we briefly discuss the condensation energy 
Uq. In the case of the pure DSC state, electrons 
near (tt, 0) contribute much to Uo and those near 
(ir/2, ir/2) contribute little. However, in the DDW sce- 
nario, there is a pre-existing Wddw in the pseudogap 
region by assumption. Upon doping, the Fermi sur- 
face is a small pocket, which has not been seen in 
experiments yet. The same d x 2_ y 2 symmetry makes 



the vicinity of (n, 0) already far below the Fermi sur- 
face. Then paring can not affect them as significantly 
as in the case of pure DSC, thus Uo is reduced. Let 
us study the contribution to Uo along the curve of 
the minimum gap form the direction (tt/2,tt/2) — > 
(-7T, 0). At beginning, it is proportional to Ansc^fej 
where (f>k = cos k x — cos k y . After passing the point 
kg where the end of the Fermi pocket lies, the depen- 
dence changes into ^(a/^ddwI 1 = <feo/0fc) 2 + A^ sc - 
W 7 ddw(1 — 0fco/0fe))- If Wddw is large, the slope becomes 
softer and an apparent kink develops. This kink can be 
testified by studying the retreat of the leading edges of 
the high resolution ARPEPS data deep in the supercon- 
ducting region relative to those in the pseudogap region. 
If the pseudogap is caused by pair fluctuations, such kink 
will not exist. 

In summary, we studied the DDW and DSC order 
parameters' dependence with temperature and thermo- 
dynamic quantities in detail by solving the mean-field 
Hamiltonian self-consistcntly. The DDW order is sup- 
pressed when temperature drops below T c in the under- 
doped region because of their competing nature. The 
disagreement with experimental results was discussed. 
Behaviors of the chemical potential, entropy and spe- 
cific heat with temperature and doping are investigated. 
The increase of chemical potential is predicted when the 
temperature increases in the pseudogap region. We also 
showed the decrease of entropy when doping decreases in 
the underdoped region. The distribution of condensation 
energy in the momentum space has a kink along the di- 
rection form (7r/2, 7r/2) to (7r,0). These features may be 
used in experiments to testify the DDW scenario. 
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